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We consider the Higgs sector of multi-Higgs-doublet models in the presence of simple symmetries 
relating the various fields. We construct basis invariant observables which may in principle be used 
to detect these symmetries for any number of doublets. A categorization of the symmetries into 
classes is required, which we perform in detail for the case of two and three Higgs doublets. 
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, Many features of the Standard Model (SM) of electroweak interactions have been accurately tested. Still, the Higgs 
• sector remains largely unknown. Indeed, even after one scalar particle is directly detected, there may well be further 
p\| ' scalars awaiting discovery (as required, for example, by supersymmetry) . It is easy to construct a iV-Higgs-doublet 
model (NHDM), but the number of parameters in the Higgs potential grows very rapidly with N. A generic two- 
Higgs-doublet model (THDM) has 14 real parameters in the Higgs potential, while the generic three-Higgs-doublet 
model (3HDM) already has 54 real parameters. 
I ' The number of parameters will be reduced if the theory has discrete (or continuous) symmetries relating the 
various Higgs fields, which we denote by Higgs Family symmetries or HF-symmetries. Besides parameter reduction, 
such symmetries may also be desirable features of a theory in order to preclude flavor changing neutral currents or to 
explain relations among different observables. This article presents some features of HF-symmetries in the NHDM. 

Two deceptively simple questions about HF-symmetries arise. First, classifying the symmetries by their impact on 
the Higgs potential, one would like to know how many distinct classes of symmetries may be implemented. Surprisingly, 
this turns out to be a rather nontrivial question. Second, a simple basis change among the various Higgs fields alters 
the Lagrangian but, obviously, not its physical consequences. Signals of HF-symmetries should be invariant under 
these transformations. 
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fS| ' The need to seek basis invariant observables in models with many Higgs was pointed out by Lavoura and Silva 
I [l[, and by Botella and Silva 0, stressing apphcations to CP violation. Refs. 0,9 indicate how to construct basis 
• invariant quantities in a systematic fashion for any model, including multi-Higgs-doublet models. Work on basis 
QQ ' invariance in the THDM was much expanded upon by Davidson and Haber by Gunion and Haber 0], and 
by Haber and O'Neil JJ. Basis invariance in the THDM was also considered in Refs. i, [m, El- In particular, 
^ Davidson and Haber Hldevelop several strategies to construct basis invariant descriptions of HF-symmetries, in the 
. ^ context of the THDM One of our aims is to extend their work into multi-Higgs systems. 

' The paper is organized as follows. In section [IT] we introduce our notation and show that a simple HF-symmetry 
H I may always be reduced to a standard diagonal form through a basis transformation. Then, we turn to the problem of 
- - ' classifying the HF-symmetries according to their action on the Higgs potential. We cover the THDM in section [1111 
and we discuss the 3HDM in section IIVI In section |V] we define a set of basis invariant observables applicable to 
any NHDM, which may in principle be used in order to identify the presence of HF-symmetries. We present our 
conclusions in section IVTl Appendix [X] includes the implications that the different classes of symmetries of the 3HDM 
have on the quadratic and quartic coupling coefficients. This provides the fingerprint database against which the basis 
invariant observables of section |V] should be compared. 
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II. THE SCALAR SECTOR OF A GENERIC N-HIGGS-DOUBLET MODEL 
A. The scalar potential and basis transformations 

In this article we follow the notation of Refs. 0, 0, H]- Let us consider a SU{2) (g) U{1) gauge theory with N 
Higgs-doublets $i, with the same hypercharge 1/2, and with vacuum expectation values (vevs) 







(1) 



The index i runs from 1 to N, and we use the standard definition for the electric charge, whereby the upper components 
of the SU{2) doublets are charged and the lower components neutral. 
The scalar potential may be written as 



where Hermiticity implies 



Y — Y* 



(2) 



(3) 



The number of independent parameters of this potential is shown in Table U 

TABLE I: Number of parameters in the Y and Z coefficients of the Higgs potential. 
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The stationarity conditions are 

[Y^j + 2Zij^ki vlvi] Vj = 
Multiplying by v* leads to 

Ytjiv*iyj) = -2Zy- fc; {v*Vj) [vlvi) 



(fori = ,N). 



(4) 
(5) 



We may rewrite the potential in terms of new fields $^ , obtained from the original ones by a simple basis transfor- 
mation 



(6) 



where U is a N x N unitary matrix. Under this unitary basis transformation, the gauge-kinetic terms remain the 



same but the coefficients Yij and Zij^ki are transformed as 



and the vevs are transformed as 



^ij.kl ^ ^ij,kl — ^im Uj^o ^mn,op ^jn ^Ipi 



(7) 
(8) 



(9) 



Thus, the basis transformations U may be utilized in order to absorb some of the parameters in Y and/or Z, meaning 
that not all parameters in Table |l] have physical significance. 
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B. Higgs Family symmetries 

Let us assume that the scalar potential in Eq. ([2]) has some explicit internal symmetry. That is, we assume that 
the coefficients of Vh stay exactly the same under a transformation 

^ $f = 5.,$,. (10) 

S' is a unitary matrix, so that the gauge-kinetic couplings arc also left invariant by this HF-symmetry. As a result of 
this symmetry 

y^, - r^^ = s.kYkiS*i, (11) 

Zij^kl — ^ij,kl ~ ^im Sko Z„in,op S Sip, (12) 

Notice that this is not the situation considered in Eqs. ©-([S]). There, the coefficients of the Lagrangian do change. 
What we said there was that, although the coefficients do change, the quantities which are physically measurable 
cannot. What we consider in Eqs. (fT0 |) -(fT2 |) is different. Here we consider the possibility that Vh has some HF- 
symmetry S which leaves the coefficients unchanged. 

We now turn to the complicated interplay between HF-symmetries and basis transformations. Let us imagine that, 
when written in the basis of fields <i>i, Vh has a symmetry S. Then we perform a basis transformation from the basis 
3>i to the basis as given by Eq. Clearly, when written in the new basis, Vh does not remain invariant under 
S. Rather, it will be invariant under 

S' = USUI. (13) 

As we change basis, the form of the potential changes in a way which may obscure the presence of a HF-symmetry. 
Eq. (jl3p means that many HF-symmetries which might look distinct on the surface, will actually imply exactly 
the same physical predictions. Any two symmetries S and 5" related by Eq. (fT3| . for some basis transformation 
U, will make the same predictions. Now S, S', and U are all matrices of the U{N) group, within which Eq. (fT5|) 
constitutes a conjugacy relation. Thus, Eq. (jl3p means that HF-symmetries associated with matrices S and 5' in 
the same conjugacy class of U{N) correspond to the same model. This result is easy to generalize because an overall 
phase transformation on U oi S has no impact on the potential Vh. This can be seen directly from Eqs. 
and Eqs. pT|) ~ (fT2)) . and is due to the fact that the Higgs potential Vh in Eq. ([2]) only depends on the Higgs fields 
through bilinear combinations. So, symmetries S and S' belonging to conjugacy classes related by a global phase 
transformation lead to the same physics. 



C. A special basis 



One can show that a N x N complex matrix S belongs to U{N) if and only if there exists a unitary matrix U 
such that S' in Eq. is diagonal, with all entries of magnitude 1 [l^. This means that, by a suitable basis 
transformation, any symmetry S may be brought to the form 



/ 



(14) 



where < 9i < 27t {i — 1 . . . N). The conjugacy classes can thus be classified by matrices of the type in Eq. 
In the basis where the symmetry is represented by Eq. (|14p . the coefficients must obey 



■'ij.kl 



ij.kl ; 



(15) 
(16) 



where there is no sum over repeated indexes. This result is obtained by substituting the special form of S in Eq. JT^ 
onto Eqs. (|li p - (fT2l) . The result in Eq. ([15]) applies not only to the matrix Y but to any matrix whose two indexes 
transform as U^ and U, as shown for Y in Eq. ([7|). In particular, in this special basis, the matrices 



.(1) 



.(2) 



^ ^ ^ik,kj 1 
k 

^ ^ ^ij,kk ) 



(17) 
(18) 
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introduced by Davidson and Haber [J] must also obey a relation like Eq. p5|) . (Although sum over repeated indexes 
is assumed unless explicitly stated, we have shown it here explicitly for clarity.) 

For symmetries corresponding to 6i ^ 9j for all i j, Eq. (fT5|) implies that the matrix Y is diagonal. If the Higgs 
potential only had quadratic terms, this information would be useless, since any hermitian matrix can be diagonalized 
by a suitable unitary basis change. Said otherwise, without quartic terms in the potential, imposing HF-symmetries 
(or not) would make no difference. Thus, any sign of HF-symmetries must necessarily involve the quartic terms. For 
the special case of the THDM, this can be seen explicitly in Eqs. (39)-(50) of Ref. [j 

But given two matrices (for example, A = Y and B = Z'-^'^), Eq. psp already gives crucial information in the case 
where 6i 7^ 9j for all i ^ j. Indeed it states that, in the special basis, A and B are simultaneously diagonal. As a 
result, their commutator vanishes; [A, B] = 0. Now, the commutator is a matrix and the null matrix is always mapped 
onto the null matrix, regardless of which basis transformation one chooses. Thus, we conclude that symmetries which 
are represented in the special basis by 9i ^ dj (for all i ^ j) will lead to the basis- invariant result [A, B] = 0. This 
can be used to define basis invariant fingerprints of HF-symmetries, explaining Eqs. (39)-(41) of Ref. [j] 

One could now ask whether all possible impositions due to HF-symmetries can be cast in the form [A, B] ^ for 
suitably chosen matrices A and B. The answer is negative, as Davidson and Haber found when trying to disentangle 
the Peccei-Quinn [14J symmetry from the usual Z2 symmetry [isj - c.f. their Eq. (46). 



D. Further simplifications due to global phase invariance 



Because the Lagrangian is invariant under global phase transformations, there are infinitely many conjugacy classes 
which imply the same physical predictions. Indeed, classes represented by the diagonal elements 



( 1 



\ 



\ 



JO2 



(19) 



that we can concentrate on symmetries of the type shown in Eq. ((T9)l . without the pre- factor 
shall classify each class of symmetries by their diagonal representative: 



for fixed values of 9j {j = 2 . . . N) lead to the same physical predictions, regardless of the value of 9i. This means 

Thenceforth, we 



(20) 



/ 1 



V 



Alternatively, we could use the e*^^ phase freedom in order to restrict our attention to symmetries in SU{N). 

A first question now arises: do classes corresponding to different values of 9j {j = 2... N) necessarily imply different 
physical predictions? The answer is negative. A second question arises: can one classify the different types of HF- 
symmetries according to their impact on the Higgs potential? The answer is affirmative, but that must be done 
separately for each value of N. We review the case of = 2 in Section UIIl and we turn to the more difficult case of 
A = 3 in Section HYl 



III. HF-SYMMETRIES IN THE THDM 



In the previous sections we learned the following. Two symmetries in the same conjugacy class yield the same 
physics. Thus, we can go into a special basis and consider a diagonal matrix with complex entries of unit magni- 
tude. In fact, there are infinitely many such diagonal matrices which yield the same physics, because global phase 
transformations have no impact on the Lagrangian. Therefore, we can concentrate on symmetries of the type 



(21) 



5 



Substituting into Eq. (|T5|) . we obtain 




We conclude that the Y matrix elements come affected by the following phase factors: 

-a 
a 



(22) 



(23) 



This table is a shorthand notation to keep track of the exponents which appear in Eq. (j22p . If S" is indeed a symmetry 
of the potential, then these phase factors must equal (mod 27r) for any nonzero value of the corresponding entry in 
the Yij matrix. If a = 0, we have the uninteresting identity transformation. We shall ignore this possibility henceforth. 
If < a < 27r, then any matrix in the problem (built from the coefficients in the scalar potential) must be diagonal. 
This leads to conditions of the type [A, B] = discussed above. Notice that this condition does not distinguish a — n, 
corresponding to the Z2 symmetry 



^1 




(24) 



-a 




—a —2a 


a 




-a 



from the transformation with a = tt/S, etc. Accordingly, Davidson and Haber were unable to find a condition to 
distinguish Z2 from Peccei-Quinn based exclusively on matrix conditions of the type [A, B] = 0. 

The Z2 symmetry will only be distinguished from the symmetries with other values of a by the quartic terms. 
Substituting Eq. (|2ip into Eq. we conclude that the Z matrix elements come affected by the following phase 

factors: 



(25) 



This is represented as a table of tables. The uppermost-leftmost table corresponds to the phases affecting Zn ^i. The 
next table along the same line corresponds to the phases affecting Z12M, f^nd so on... If S is indeed a symmetry 
of the potential, then these phase factors must equal (mod 27r) for any nonzero value in the corresponding entry 
of the Zij ^.1 tensor. Unlike what happened for the quadratic terms (and, in general, for any matrix built out of 
quadratic and/or quartic terms) we see that there is a distinction between two cases, according to whether 2a = 2tt 
or 2a ^ 0, 27r. If a = tt, then the terms 2^12,12 and Z21.21 (which are related to a parameter denoted by A5 in usual 
presentations of the THDM) may be different from zero. In contrast, A5 must vanish for symmetries with a 7^ 0,7r. 
So, the quartic terms do distinguish Si in Eq. (j24ll from 



a 




-a 


2a a 




a 




(26) 



But they do not distinguish among the symmetries 





(27) 



These symmetries are actually quite curious. Suppose we impose the symmetry S'2/3 on the Lagrangian. Clearly, 
applying the symmetry again must also leave the Lagrangian invariant. As a result, the Lagrangian is invariant under 



'2/3: 



02 

•^2/3' 



and Sl/^ 



= 1, which form a closed group. The Lagrangian is always invariant under a group; if it is 
invariant under symmetries Sa and Sh, it is obviously also invariant under SaSt- And if we choose a/ir irrational, the 
group 5*2, iS*!, S'l, . . . will even have an infinite number of elements. 

But there is a further important point. We have just shown that if the potential is invariant with respect to a 
symmetry S2 for some value of a ^ 0, tt, then it will necessarily be invariant with respect to a symmetry S2 with 
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any other value of a. That is, we have imposed a discrete symmetry but the resulting potential is invariant with 
respect to a continuous symmetry - the Peccei-Quinn symmetry (14l] . This is an important point because continuous 
symmetries, if broken, imply the presence of massless Goldstone bosons. Suppose we build a NHDM with an innocent- 
looking discrete symmetry. It may happen that imposing this symmetry has the same effect on the potential as a 
global symmetry and, thus, the possibility exists for undesired massless scalars. We have just seen one such example. 
We impose the discrete symmetry 5*2/3 (the corresponding group of symmetries, to be precise) only to find that the 
resulting potential is invariant under the continuous Peccei-Quinn symmetry. 

Notice that the analysis of the quartic terms is sufficient to isolate all cases of interest. Indeed, the uppermost- 
leftmost 2x2 block of Eq. ((25)) coincides with Eq. ((23|) . A similar situation occurs for any other value of N. 

In conclusion, as far as simple HF-symmetries are concerned, we have only three possibilities. Either we have 
the most general Lagrangian, or we have the Z2 symmetry, or we have the PQ symmetry. This exhausts all simple 
symmetries. We are not considering here CP-type symmetries, and we comment briefly on multiple symmetries in 
section IIII Al 



A. Multiple symmetries 

In this article we concentrate on what we call simple symmetries. By this we mean the following: we choose some 
symmetry S and we impose only that symmetry on the Higgs potential. 

We recall two points. First, under a basis change, that symmetry will look different. For example, if we impose 
the symmetry $1 $1, $2 ~> — on some basis, then that symmetry will turn into ^[ ^ $2 if we change into the 
basis = ($1 + $2)/v^, $2 = (*i - *2)/\/2. Indeed, 

is the corresponding Eq. (fT3|) 

Second, it may be that imposing that symmetry alone will yield a potential with a larger symmetry. We are not 
referring only to the obvious possibility that the potential becomes automatically invariant to the group of all powers 
of S (as exemplified above in connections with >5'2/3). It may be that the potential becomes automatically invariant 
under other symmetries, such as continuous symmetries, as indeed happens in the THDM. 

All the possibilities discussed thus far fall under what we call simple symmetries. But we may have more complicated 
situations. We may impose simultaneously two symmetries. For example, we may ask that the potential be invariant 
both under 5*1 in Eq. (|24p and D in Eq. (j28p . in the same basis. Notice that now it is irrelevant that 5*1 and D are 
in the same conjugacy class when considered individually. We are imposing both in the same basis. Bringing one 
to diagonal form will make the other off-diagonal, and vice-versa. In this case, the potential becomes automatically 
invariant under the group of four symmetries 5*1, D, SiD, and 1. This case is considere d by Davidson and Haber 
Q after their Eq. (37). Other multiple symmetries of the THDM were studied by Ivanov [I4I. The general analysis 
of such cases in the NHDM is much more difficult and it is not considered in this article. Indeed, what we dubbed 
simple symmetries will prove surprisingly demanding, even for = 3. 



IV. HF-SYMMETRIES IN THE 3HDM 



The analogues of Eqs. fUl), (ESI), and ([251), for iV = 3 are 




-a -(3 
a a — P 
(3 13-a 



(30) 
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and 



-a -(3 
a a — P 
(3 P-a 

a a — (3 
2a a 2a — P 

a + P P a 

P P-a 
a + P P a 
2/3 2P-a P 



—a —2a —a — P 
-a -P 
P — a P — 2a —a 



-a -P 
a a — P 
P P-a 



P — a P — 2a —a 

P P-a 
2/3 -a 2/3 -2a P-a 



-P -a-P -2/3 
a-P -P a - 2/3 
-a -P 



a-P 



a - 


-2/3 


2a 


-2/3 
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-P 



-a -/3 
a a-P 
P P-a 



(31) 



respectively. 

Let us first look at the impact of the symmetries on the quadratic terms in Eq. (j30p . One interesting situation 
is a = 0, /3 ^ 0. In this situation <i>i and $2 have the same transformation, while <i>3 transforms differently. One 
could think of a 7^ 0, /3 = as a different situation, but it is not. It is the same as the previous situation, with the 
interchange of fields 2 and 3. Since such a field permutation corresponds to a basis change (achievable through some 
unitary matrix U), the two situations correspond to exactly the same symmetry viewed in different basis, and lead 
to the same physics. Now we consider a ^ 0, /3 7^ 0. The (2, 3) and (3, 2) entries in Eq. ([50)1 show that we must 
distinguish a = P from a ^ p. But a — P corresponds to the symmetry 



1 



v 



(32) 



1 



which, aside from the irrelevant overall phase, is just a 1 <-> 3 permutation of the situation already considered, 
conclusion, the quadratic terms in the Higgs potential distinguish among the following two types of symmetries: 



In 



/ 



V 



1 



/ 



V 



1 



(33) 



pi/3 



As happened for N = 2, looking at the quartic terms will further open up these classes of symmetries. Said 
otherwise, two symmetries may have the same impact on the quadratic terms but different impact on the quartic 
terms. The zeros in Eq. (jSlJ correspond to the entries of the Z tensor which are real. In order to see how the 
symmetries affect the quartic terms we start by collecting all distinct combinations in Eq. (|3ip : a, 2a, /3, 2/3, a + P, 
a-P, 2a — /3, a — 2/3, and 2a — 2/3. Each may be equal to (mod 2tt) or not. We study all possible combinations, 
making sure that each new class of symmetries found does not correspond to a mere basis transformation of a class 
considered previously. Proceeding in this fashion, we find that (as far as simple HF-symmetries are concerned, and 
aside from the most general Lagrangian) there are seven types of symmetries having distinct impacts on the Higgs 
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potential: 



^1 



-1 



^2 



o2i7r/3 



-2i7r/3 



/ 1 



^4 = 



^6 = -1 



ia=^0,-!T/2,2TT/3,Tr) 



S7 = 



(a/0,7i-|/3^±a,0,7r) 



(34) 



Comparing with Eq. (|33|) . we see that 5*1 and S2 have the same impact on the quadratic terms. This is different from 
the impact of 83-87 on the quadratic terms. The impact of the symmetries 8i~8i on the coefficients of the Higgs 
potential is presented in appendix [Xl 

Recall that the potential is always invariant under a group of symmetries. For example, imposing S3, the potential is 
automaticaly invariant under 83, S3, and 83 — 1. The symmetries 82, 85, and 5*6 are special in this respect. Imagine 
that we impose 82 from some value of a 0, tt. Then, the potential is automaticaly invariant with respect to 5*2 with 
any other value for a. That is, we wish to impose a discrete symmetry, but the resulting potential turns out to be 
invariant under a J7(l) continuous symmetry. The case for 87 is even worse. Imposing 87 for some chosen numerical 
values for a (7^ 0, tt) and /3 zta, 0, tt) will automaticaly generate a potential invariant under all symmetries 6*7 for 
any values of a and /?. The resulting potential will be invariant under a U{1) (g) U{1) continuous symmetry. 



V. BASIS INVARIANT DESCRIPTIONS OF THE SYMMETRY CLASSES IN THE NHDM 

Let us look back at Eq. (|33p. We recall that the analysis of the quadratic terms applies equally well to any matrix 
(even if built out of the quartic terms, as Z^^^ and Z'^'). Therefore, a (basis invariant) commutator type condition 
will distinguish one symmetry of the first type in Eq. ([55)1 {81-82, where [A, B] 7^ 0) from one of the second type 
{83-87, where -B] = 0). But it will not distinguish among two symmetries of the same type in Eq. ([55]) . For 
example, it will not distinguish 83 from 87, even though they have a different impact on the quartic terms. As a 
result, commutator conditions, which were so central to Davidson and Haber's study of the THDM [4|, have a very 
limited use for TV > 3. 

Basis invariant fingerprints of the HF-symmetries may be found by combining eigenvectors of the matrix of quadratic 
couplings Y with the quartic couplings Z. In the context of the THDM this seemed a curiosity and was left by Davidson 
and Haber Q to the end of their appendix B; Eqs. (B17)-(B22). Inspired by this remark, we developed a technique 
which, when suitably extended and interpreted, will become central to our definition of the basis invariant observables 
identifying HF-symmetries. 

Consider a general 3HDM. We define the three eigenvectors of the matrix Y by y^, y^, and y^. If the matrix has two 
or three degenerate eigenvalues, we are free to choose any orthonormal basis in the degenerate space. The components 
of the eigenvector are denoted by yf . Under a basis change U, the components of the eigenvectors change as 

vT - u., vr (35) 

Combining this with Eq. ([8]) we see that the quantities 

i-P'-<'^z,,^u{y?Y{yl){yir{yf) (36) 
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are basis invariant for any values of a, /9, 7, and S between 1 and 3. Therefore, we can evaluate them in any basis. In 
particular, in a basis where Y is diagonal, we may choose the eigenvectors as 



That is, y" — Sai, where Sai is the Kronecker symbol. In the basis of Eq. ([57)) . the quantities in Eq. become 



N'^{N'^ + l)/2 of these are independent. 

Before we proceed, we must point out a subtlety concerning Eqs. (|36p - (|38p . Suppose that we have a Y matrix 
whose eigenvalues are not degenerate, and that we find its three eigenvectors. Now we have a problem in attributing 
to them the labels 1, 2, and 3. There are 6 possibilities, which differ by permutations. What one author chooses as 
/i,i,i,2 ^^-y |-,g -^jjat another author chooses as 72,2,2,1 Qj-^^-g ^jjjg choice is made, then the quantity is basis invariant. 
But the choices of different authors may differ by permutations connected with their specific choices for the ordering 
of the eigenvectors. Eq. ([55]) identifies one possibility; the other five possibilities differ by permutations in the choice 
of basis eigenvectors in Eq. ([37]) Thus, permutations must be considered when using appendix]^ as explained below. 

We now show how the quantities in Eq. p6p can be used to identify the various discrete symmetries. Let us assume 
that the Higgs potential is invariant under some symmetry S-^-S^; for example S'3. The potential may be originally 
written in a basis where this symmetry is not diagonal. But that is irrelevant; we may always consider what happens 
in a basis where the symmetry has the form in Eq. (|34p . In this basis the Z coefficients have the structure (of zero 
and non-zero entries) presented in appendix VK\ for the S'3 symmetry. Also in this basis, the Y matrix is diagonal. 
Therefore, its eigenvectors are given by Eq. ([57]) or some permutation thereof. We conclude from Eq. ([55)) that the 
observables in Eq. p6p must fall into the pattern shown in appendix \K\ which corresponds to S'3, or some permutation 
thereof. 

The discussion in the previous paragraph invoked a special basis, only to show that the presence of a symmetry 
S3-S7 will force the observables of Eq. p6p to fall onto the corresponding pattern shown in appendix |^ (or some 
permutation thereof). But the invariant need not be calculated in this basis. Because it is a basis invariant, it can 
be calculated in any basis whatsoever; the result must be the same. So, the algorithm to identify the presence of a 
symmetry is straightforward: 

• We start with the potential in some original basis. The potential has the symmetry Si in that basis (in general 
Si will not have the simple diagonal form when written in that original basis). 

• We find the eigenvectors of Y (which, in general, will also not be diagonal in the original basis). 

• We combine the Y eigenvectors with Z to calculate the basis invariant I"f^''>^ observables in Eq. (|36p. 

• We check whether the resulting pattern matches the patterns in appendix \X\ (or some permutation thereof). 

This procedure identifies which symmetry we have, even when the potential is written in an original basis where Si 
has a very obscure form. 

We have postponed the proof that this procedure also works for Si {i — 1,2) because the basis where this symmetry 
is diagonal does not guarantee that Y is diagonal. Indeed, in a basis where Si (say) is diagonal, Y is block diagonal, 
c.f. appendix [X] In order to make Y diagonal and guarantee that its eigenvectors can be cast in the form of Eq. ([ST]) 
(aside from permutations), we need a further diagonalization of the uppermost-leftmost 2x2 block of Y . But because 
the two first eigenvalues of Si and S2 are degenerate, a unitary 2x2 rotation on the uppermost-leftmost block has 
no effect on the form of the symmetry. This shows that a basis may be found where S1-S2 have the form in Eq. p4p 
and Y is diagonal, completing our proof. 

Clearly, Eqs. (|36p and ([55]) hold for any value of iV. As a result, we have succeeded in defining basis invariant 
quantities which can in principle be utilized in order to identify any HF-symmetry in NHDM, for any value of N. 
But in order to perform this identification in practice we need to have a set of textures to compare with, as we have 
done in appendix [X] for = 3. The problem of categorizing the different classes of HF-symmetries which may affect 
the Higgs potential (and, thus, the corresponding textures) becomes demanding as A^ increases. For example, in a 
cursory analysis of A^ 4 we have identified at least 15 distinct classes of symmetries. 




(37) 
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One final remark concerns the possiblitity that the matrix Y has degenerate eigenvalues. In this case we must define 
new I'^-'l^'i'^ parameters invoking the eigenvectors of Z^^^ (or Z'^^^) rather than the eigenvectors of Y. Some regions 
of parameter space may require special care. These types of questions were already present in the various methods 
proposed for the THDM ^. 

An apt analogy to our procedure is the following. We wish to identify a symmetry. The I°'f^-i^ in Eq. ((36)) provide 
us with a (basis independent) fingerprint of the symmetry. But we can only use this information in order to identify 
the symmetry, if we have a database with all the distinct fingerprints which may show up, one for each symmetry class. 
This is what we provide explicitly in appendix for iV = 3. Anyone interested may construct a similar database for 
iV > 4. 



VI. CONCLUSIONS 

We have constructed a set of basis invariant quantities which may in principle be used to identify the presence of 
HF-symmetries in a NHDM, regardless of the value of N . HF-symmetries can be classified according to their impact 
on the Higgs potential. Surprisingly, this classification is already involved for iV = 3. We have discussed the cases of 
N = 2 and = 3 in detail showing how to combine the J"^'T* observables with the classification scheme in order to 
identify any HF-symmetry, regardless of the basis in which the Higgs potential may be originally written in. Our basis 
invariants may be applied to any other value of A^, by constructing the database of the classes of symmetries 

possible for that value of A^. 

This classification is also important because sometimes one imposes a discrete symmetry only to find that the 
potential becomes automatically invariant under a much larger class of symmetries. These may even be continuous, 
implying the danger that Goldstone bosons might appear. We provide explicit examples of this problem. This 
may even be relevant for studies of the fermion sector. For example, Grimus et al. discuss very general symmetry 
realizations of texture zeros in the fermion sector with the help of scalar fields onto which certain discrete symmetries 
are imposed [Tg} . When using such techniques, one must inspect also the Higgs potential in some detail, including 
the symmetry breaking, lest there be undcsired massless scalars. 
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APPENDIX A: COUPLING STRUCTURES FOR THE DIFFERENT CLASSES OF SYMMETRIES 



In this appendix we discuss the impact that the seven classes of symmetries identified in the 3HDM have on the 
couphng constants in the scalar potential. We show the result in the special basis in which the symmetry has one of 
the diagonal forms in Eq. p4p . In section |V] we show how to turn this information into a basis invariant fingerprint 
for the discrete symmetries. 

The quadratic couplings distinguish three cases: i) the most general potential, where all entries of Yij may be 
nonzero; ii) the potential with one of the symmetries S'i-5'2, where the matrix Yij is block diagonal and the uppermost- 
leftmost 2x2 block is left unconstrained; and iii) the potential with one of the symmetries 83-87, where the matrix 
Yij is diagonal. The first case corresponds to 3 real and 3 complex parameters (for a sum of 9 real variables); the 
second case corresponds to 3 real and 1 complex parameters (5 real variables); and the third case corresponds to 3 
real parameters. 

To see the impact on the quartic potential, we organize the Zij^ki tensor into a matrix of matrices. The uppermost- 
leftmost matrix corresponds to the phases affecting Zn^ki- The next matrix along the same line corresponds to the 
phases affecting Zi2^ki, and so on... We use the following notation for the various entries 
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(Al) 



where (i = 1 . . .9) are real and q (i = \ . . . 18) are complex. In the basis of Eq. ([M)) . the symmetries 5'i-S'7 set 
different combinations of Ci to zero but leave the real coefficients unconstrained. 

The most general 3HDM has 9 real and 18 complex quartic couplings for a total of 45 real variables. Combining with 
the quadratic parameters, we have 12 real and 21 complex parameters, for a total of 54 real variables (33 magnitudes 
and 21 phases). However, not all the variables have physical significance due to the possibility of changing basis 
through any 3x3 unitary matrix (which can be parametrized with 3 magnitudes, 5 relative phases, and 1 global 
phase). Thus, the most general 3HDM has 30 magnitudes and 16 phases with physical significance. 

If the potential obeys the symmetry 8\ , in the basis of Eq. ([M)) the quartic couplings have the following structure, 
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This corresponds to 9 (12) real and 8 (9) complex parameters in the quartic couplings (in the scalar potential). 
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If the potential obeys the symmetry S2, in the basis of Eq. ((34)) the quartic couphngs have the following structure, 
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This corresponds to 9 (12) real and 5 (6) complex parameters in the quartic couplings (in the scalar potential). 
If the potential obeys the symmetry 53, in the basis of Eq. p4p the quartic couplings have the following structure, 



ri 
















C4 







C4 










r4 







r7 
















Cll 










rs _ 







C12 




















r7 







r4 
















Cll 










r* 
C12 







r2 







r* 
C12 










r* 

C4 
















re _ 







rg 













rs 







C12 







rs 










c| 
















rg 







re 










Cll 







r* 
Cll 
















r3 _ 





(A4) 



This corresponds to 9 (12) real and 3 (3) complex parameters in the quartic couplings (in the scalar potential). 
If the potential obeys the symmetry ^4, in the basis of Eq. (|34p the quartic couplings have the following structure. 
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This corresponds to 9 (12) real and 2 (2) complex parameters in the quartic couplings (in the scalar potential). 
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If the potential obeys the symmetry S5, in the basis of Eq. ((34)l the quartic couphngs have the following structure, 
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This corresponds to 9 (12) real and 1 (1) complex parameters in the quartic couplings (in the scalar potential). The 
single complex parameter appears in four entries. 

If the potential obeys the symmetry Se, in the basis of Eq. (|34p the quartic couplings have the following structure, 
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This corresponds to 9 (12) real and 1 (1) complex parameters in the quartic couplings (in the scalar potential). Unlike 
what happened for 5*5, here the single complex parameter appears in two entries. 

If the potential obeys the symmetry Sj, in the basis of Eq. ([34| the quartic couplings have the following structure. 
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This corresponds to 9 (12) real parameters in the quartic couplings (in the scalar potential). All complex parameters 
in the scalar potential vanish. 



